integrate (alpha+beta*cosh(B*t)^2)/(gamma+delta*cosh(B*t)^2) dt